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MHD SIMULATIONS OF BONDI ACCRETION TO A STAR IN THE "PROPELLER" REGIME 
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ABSTRACT 

This work investigates Bondi accretion to a rotating magnetized star in the "propeller" regime using 
axisymmetric resistive, magnetohydrodynamic simulations. In this regime accreting matter tends to be 
expelled from the equatorial region of the magnetosphere where the centrifugal force on matter rotating 
with the star exceeds the gravitational force. The regime is predicted to occur if the magnetosphcric 
radius larger than the corotation radius and less than the light cylinder radius. The simulations show 
that accreting matters is expelled from the equatorial region of the magnetosphere and that it moves 
away from the star in a supersonic, disk-shaped outflow. At larger radial distances the outflow slows 
down and becomes subsonic. The equatorial matter outflow is initially driven by the centrifugal force, 
but at larger distances the pressure gradient becomes significant. We find that the star is spun-down 
mainly by the magnetic torques at its surface with the rate of loss of angular momentum L proportional 
to , where f2* is the star's rotation rate and /i is its magnetic moment. Further, we find that 

L is approximately independent of the magnetic diffusivity of the plasma r\ m . The fraction of the Bondi 
accretion rate which accretes to the surface of the star is found to be cx fij 1 - // -1 - 7 ^ 4 . Predictions of 
this work are important for the observability of isolated old neutron stars and for wind fed pulsars in 
X-ray binaries. 

Subject headings: accretion, dipole — plasmas — magnetic fields — stars: magnetic fields — X-rays: 
stars 



1. INTRODUCTION 

Rotating magnetized neutron stars pass through differ- 
ent stages in their evolution (e.g., Shapiro & Teukolsky 
1983; Lipunov 1992). Initially, a rapidly rotating (P <J Is) 
magnetized neutron star is expected to be active as a ra- 
diopulsar. The star spins down owing to the wind of mag- 
netic field and relativistic particles from the region of the 
light cylinder (Goldreich and Julian 1969). However, af- 
ter the neutron star spins-down sufficiently, the light cylin- 
der radius becomes larger than magnetospheric radius r m 
where the ram pressure of external matter equals the mag- 
netic pressure in the neutron star's dipole field. The rel- 
ativistic wind is then suppressed by the inflowing matter 
(Shvartsman 1970). The external matter may come from 
the wind from a binary companion or from the interstellar 
medium for an isolated neutron star. The centrifugal force 
in the equatorial region at r m is much larger than grav- 
itational force if r m is much larger than the corotation 
radius r cor . In this case the incoming matter tends to be 
flung away from the neutron star by its rotating magnetic 
field. This is the so called "propeller" stage of evolution 
(Davidson & Ostriker 1973; Illarionov & Sunyaev 1975). 

The "propeller" stage of evolution, though important, 
is still not well-understood theoretically. Different studies 
have found different dependences for the spin-down rate 
of the star (Illarionov & Sunyaev 1975; Davis, Fabian, & 



Pringle 1979; Davis & Pringle 1981; Wang & Robertson 
1985; Lipunov 1992). Observational signs of the propeller 
stage have been discussed by number of authors (e.g., 
Stella, White, & Rosner 1986; Treves & Colpi 1991; Cui 
1997; Treves et al. 2000). 

MHD simulations of disk accretion to a rotating star 
in the propeller regime were done by Wang and Robert- 
son (1985). However, authors considered only equatorial 
plane and concentrated on investigation of instabilities at 
the boundary between the magnetosphere and surround- 
ing medium. Thus they could not investigate accretion 
along the magnetic poles of the star. An analytical model 
of disk accretion in the propeller regime was developed by 
Lovelace, Romanova, and Bisnovatyi-Kogan (1999). Disk 
accretion at the stage of weak propeller (r TO ~ r cor ) inves- 
tigated numerically by Romanova et al. (2002). 

The mentioned studies obtained possible trends of the 
propeller stage of evolution. However, two and three di- 
mensional MHD simulations are needed to obtain more 
definite answers to the important physical questions. The 
questions which need to be answered are: (1) What are 
the physical conditions of the matter flow around the star 
in the propeller regime of accretion? (2) What is the spin- 
down rate of the star, and how does it depend on the 
star's magnetic moment and rotation rate and on the in- 
flow rate of the external matter? (3) What is the accretion 
rate to the surface of the star and how does it depend on 
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Fig. 1. — Geometry of the MHD simulation region, where Mb is the Bondi accretion rate, fi and are the magnetic 
moment and angular velocity of the star, R* is the radius of the star, (R m ax, Z max ) are the limits of the computational 
region. In the described calculations i?» <C R ma x- 



the star's rotation rate and magnetic moment? (4) What 
are the possible observational consequences of this stage 
of evolution? 

This paper discusses results of axisymmetric, two- 
dimensional, resistive MHD simulations of accretion a ro- 
tating magnetized star in the "propeller" regime. We treat 
the case when matter accretes spherically with the Bondi 
accretion rate. Bondi accretion to a non-rotating and a 
slowly rotating star was investigated by Toropin ct al. 
(1999; hereafter T99). It was shown that the magnetized 
star accretes matter at a rate less than the Bondi rate. 
Toropina ct al. (2002; hereafter T02) confirmed this re- 
sult for stronger magnetic fields. In this paper we consider 
initial conditions similar to those in T02, but investigate 
the case of rapidly rotating stars. Section 2 gives a rough 
physical treatment of the propeller regime. Section 3 de- 
scribes the numerical model and computations. Section 4 
describes the main results from our simulations, and Sec- 
tion 5 gives a numerical application of our results. Section 

6 gives the conclusions from this work. 

2. PHYSICS OF THE "PROPELLER" REGIME 

In the propeller regime the magnctosphcric radius r m is 
larger than the corotation radius r cor = (CM/f^ 2 ) 1 / 3 s=a 

0.78 x 10 9 P^ 3 cm and incoming matter in the equato- 
rial plane is flung away from the star. Here, we assume 
M = 1.4M , and P w = P/10 s with P the rotation pe- 
riod of the star. The magnetospheric radius is determined 
by a balance of the ram pressure of the inflowing matter 
Pm + Pm^m against the magnetic pressure B 2 /87r of the 
neutron star's field B m = p/r"^, where the m— subscripts 
indicate that the quantity is evaluated at a distance r m , 
and /i is the star's magnetic moment. 

Consider first the case of a non-rotating star. The Bondi 
inflow velocity at r m is supersonic for specific heat ratios 

7 < 5/3. Thus p m = M/(ATrrU 2 VGM), where M is the 
accretion rate. For this case the ram pressure is approxi- 
mately p m GM/r m so that 



pm(tt*r m ) 2 w ^ 2 /(87rr^), which gives 
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(Davidson & Ostriker 1973; Shapiro & Teukolsky 1983). 
Here, fi 30 = fi/lO 30 Gem 3 , and M_ 17 = M/T0- 17 M©/yr. 

For a rapidly rotating star, r m » r cor , the dominant 
velocity is due to the rotation of the magnetosphere so 
that w 2 m — (n*r m ) 2 . The flow velocity is much larger 
than the sound speed so that this is a "supersonic pro- 
peller" (Davies et al. 1979). Neglecting the complicated 
two-dimensional nature of the plasma flow, we estimate 
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(Wang & Robertson 1985; Lovelace, Romanova, & 
Bisnovatyi-Kogan 2002). In this limit, there is to a first 
approximation no accretion to the star; all of the incoming 
matter is flung away from the star in the equatorial plane 
by the rotating magnetic field. 

Equations (1) and (2) give the same values for SI2 = 
(2M/^ 2 ) 3 / 7 (C7M) 5 / 7 which corresponds to a period P 2 ~ 
7100(^ 2 /M_i 7 ) 3 / 7 s. The propeller regime requires peri- 
ods P < ?2- At the other extreme, the condition that r m o 
be less than the light cylinder radius gives the condition 
P > P 1 « 0.77(/4/M-i7) 2/9 s . 

The spin-down rate of the star can be estimated as 
IdQ^/dt 47rr 2 ^(r m i3 2 ^/87r), which is the magnetic torque 
at the radius r m . This assumes that B^, B m at r m . Thus 



P 



14 /13Ql 2 /13^6/13 



(3) 



dt 2rl^ 2 7 / 13 (GM) 3 / 13 ' 

This gives a spin-down time scale 

O* 2.2 x 10 7 yr 

T ~ W\ ~ , 14/13 A V6/13 p l/13 ' 

p 30 m_ 17 r w 

where the neutron star's moment of inertia / is assumed 
to be 10 45 gem 2 . The angular momentum lost by the star 
goes into the outflow of material mainly in the equatorial 
plane. The rotational power lost by the star IQ*d£l*/dt 
goes into the kinetic (thermal plus flow) energy of the 
equatorial outflow. 

3. NUMERICAL MODEL 

We simulate the plasma flow in the propeller regime us- 
ing an axisymmetric, resistive MHD code. The code incor- 
porates the methods of local iterations and flux-corrected- 
transport (Zhukov, Zabrodin, & Feodoritova 1993). The 
code is described in our earlier investigation of Bondi ac- 
cretion to a non-rotating magnetized star (T99). The 
equations for resistive MHD are 



! + V.(pv)=0 



(4) 



<9v 1 
Pg7+p(v'V)v=-Vp+-JxB|F s , (5) 



f = Vx(vxB) + ^-V 2 B, 

+ V • {pew) = -pV-v + — 
ot a 



(6) 
(7) 



We assume axisymmctry (d/d<p — 0), but calculate all 
three components of v and B. The equation of state is 
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taken to be that for an ideal gas, p — (7 — l)/?£, with 
specific heat ratio 7 = 7/5. A value of 7 less than 5/3 
is expected in the case of an ionized gas because of the 
influence of the electron heat conduction. The equations 
incorporate Ohm's law J = cr(E + v x B/c), where a is the 
electrical conductivity. The associated magnetic diffusiv- 
ity, rj m = c 2 /(47rcr), is considered to be a constant within 
the computational region. In equation (5) the gravita- 
tional force, F 9 = —GMpH/R 3 , is due to the central star. 

The star rotates with angular velocity $1* = f2* z. It 
is useful to introduce the dimcnsionlcss quantity lu* = 
n*/n K * < 1, where n K „ = {GM/Rl) 1 ' 2 is the Keple- 
rian angular velocity at the stellar radius i?* ■ Simulations 
have been done for different values of uj* = — 0.7. The in- 
trinsic magnetic field of the star is taken to be an aligned 
dipole, B = [3R(/x-R) - R 2 fi]/R 5 with /j, = p, z and 
vector-potential A = fi x R/i? 3 . 

We use a cylindrical, inertial coordinate system (r, </>, z) 
with the z— axis parallel to the star's dipole moment pi 
and rotation axis f2. The equatorial plane is treated as 
symmetry plane. The vector potential A is calculated so 
that V • B = at all times. The full set of equations is 
given in T99. 

To convert our equations to dimensionless form we mea- 
sure length in units of the Bondi radius Rb = GM/c 2 ^, 
with Coo the sound speed at infinity, density in units of 
the density at infinity poo, and magnetic field strength in 
units of Bq which is the field at the pole of the numeri- 
cal star (r = 0, z = Z*). Pressure is measured in units 
of5 2 /87r. The magnetic moment is measured in units of 
po = B Rg/2. We measure velocity in units of the Alfven 
speed vq = Bo/^/4np^. However, in the plots we give the 
velocities in units of Coo- 

After putting equations (4)- (7) in dimensionless form, 
one finds three dimensionless parameters, two of which 
are 

^-Bp 1]m =R B V = S- (8) 

Here, fj m is the dimensionless magnetic diffusivity which 
may be interpreted as the reciprocal of the Lundquist num- 
ber S; and poo = pc^ll 1S the pressure at infinity. The 
third parameter is g = GM,/(L Vg) = 7/3/2 is of the 
order of (3. We also use a parameter bo to control the 
field strength of the star. It is defined by the equation 
A = fooAo, where bo has been varied in the range 1 — 25. 

Simulations were done in a cylindrical "can" (0 < z < 
Zmax, < r < Rmax)- The size of the region was taken 
to be R max = Z max = R s /V2 = R B /hV2 « 0.141i2 B , 
which is less then the sonic radius of the Bondi flow 
R s = [(5 - 37)/ '4}R B = 0-2R B (7 = 7/5). Thus matter 
inflows supcrsonically to the computational region. The 
inflow rate is taken to be the Bondi (1952) accretion rate, 

M B = 47rA(GM*) 2 poo/Coo, where A = 0.625 for 7 = 7/5. 
The incoming matter is assumed to be unmagnetized. A 
uniform (r, z) grid with Nr x Nz cells was used. For the 
results discussed here Nr x Nz = 513 x 513. 

The "numerical star", is taken to be a cylinder of ra- 
dius i?* and half axial length Z* with i?„ = Z % << 
Rmax, Z max . In the results presented here, R* = Z. t . = 
0.0044i?B- Clearly, the numerical star is very large: the 
value of R*/Rb, although much less than unity, is much 
larger than the astrophysical ratio for a neutron star (~ 



10~ 4 ). The grid size, DR = R m ax/N R = 2.76 x 10~ 4 i? B , is 
about 16 times smaller than the star's radius. The vector- 
potential A of the intrinsic dipole field of the star is deter- 
mined on the surface of the numerical star. The value of 
the vector potential on this surface is fixed corresponding 
to the numerical star being a perfect conductor. Equiva- 
lently, the component of the B field normal to the surface 
is fixed, but the other two components vary. 

The density and internal energy at the inner bound- 
ary (the numerical star) are set equal to small numerical 
values. Thus, incoming matter with higher values of these 
parameters, is absorbed by the star. This treatment of the 
star is analogous to that used by Ruffcrt (1994). Typically, 
the density of incoming matter is ^> 10 3 times larger than 
density inside the star p = 1. For a purely hydrodynamic 
spherical accretion flow, we verified that this accretor ab- 
sorbs all incoming matter at the Bondi rate. 

4. MATTER FLOW IN THE PROPELLER REGIME 

Here, we first discuss simulations for the case = 0.5 
P = 10~ 7 , magnetic diffusivity j) m = 10~ 5 (equation 8), 
and for field strength b = 10. At the end of this section 
we discuss the dependence of the results onw„ p (and b ), 
and r\ m . 

Time is measured in units of the free-fall time from the 
top or the side of the simulation region, tff = Z max /vff, 
with Vff = \J2GM jZ max . For the rotating star, time 
is also given in terms of rotation periods of the star, 
P* = 2tt/0». 

For the mentioned values of /?, bo, and r] m , the mag- 
netospheric radius r m o and corotation radius r cor = 
(GM /n 2 ) 1 / 3, are equal for w 0.16. For smaller angu- 
lar velocities, < 0.16, the matter flow around the star 
is close to that in the non- rotating case. For > 0.16 
the flow exhibits the features expected in the "propeller" 
regime. 

Figures 2 and 3 show the general nature of the flow in 
the propeller regime. Two distinct regions separated by a 
shock wave are observed: One is the external region where 
matter inflows with the Bondi rate and the density and 
velocity agree well with the Bondi (1952) solution. The 
second is the internal region, where the flow is strongly 
influenced by the stellar magnetic field and rotation. The 
shock wave, which divides these regions, propagates out- 
ward as in the non-rotating case (T99) . For a rotating star 
in the propeller regime the shock wave has the shape of an 
ellipsoid flattened along the rotation axis of the star. 

The region of the flow well within the shock wave is ap- 
proximately time-independent. The accretion rate to the 
star becomes constant after about 1 — 2 rotation periods 
of the star. 

A new regime of matter flow forms inside the expanding 
shock wave. The rapidly rotating magnetosphere expels 
matter outward in the equatorial region. This matter flows 
radially outward forming a low-density rotating torus. The 
outflowing matter is decelerated when it reaches the shock 
wave. There, the flow changes direction and moves to- 
wards the rotation axis of the star. However, only a small 
fraction of this matter accretes to the surface of the star. 
Most of the matter is expelled again in the equatorial di- 
rection by the rotating field of the star. Thus, most of 
the matter circulates inside this inner region driven by the 
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Fig. 2. — Matter flow in the "propeller" regime for a star rotating at f2* = 0.5£1k* after 6.9 rotation periods of the star. 
The axes are measured in units of the star's radius. The background represents the density and the length of the arrows 
is proportional to the poloidal velocity. The thin solid lines are magnetic field lines. 

Fig. 3.— Enlarged view of Figure 2. The bold line represents the Alfven surface. Dotted line shows sonic surface. The 
axes are measured in units of stellar radii i?» . 

Fig. 4.— Same as on Figure 3, but the background represents angular velocity SI — i^(r, z)/r. The axes are measured in 
units of stellar radii R*. 



rapidly rotating magnetosphere. Large-scale vortices form 
above and below the equatorial torus (see Figures 2 and 
3). In three dimensions vortices of smaller scale may also 
form. Thus, at the propeller stage the significant part of 
the rotational energy of the star may go into the directed 
and thermal energy of the expelled plasma. 

The solid line in Figure 3 shows the Alfven surface, 
where the magnetic energy-density equals the thermal plus 
kinetic energy-density, B 2 /87r = e + pv 2 /2. In the equa- 
torial plane, the Alfven radius is ta ~ 4.5i?». For r < ta, 
the magnetic pressure dominates, and the magnetic field 
is approximately that of a dipole. At larger distances the 
field is stretched by the outflowing plasma. Matter inside 
the magnetopause (the region of closed magnetic field lines 
inside the Alfven surface) rotates with the angular velocity 
of the star (see Figure 4) . For r > ta, matter continues to 
rotate in the equatorial region but the angular velocity de- 
creases with r. Figure 5a shows the radial variation of the 
velocities in the equatorial plane. Matter is accelerated by 
the rotating magnetosphere for r > 4i?» up to r ~ 10i?», 
but at larger distances the radial velocity decreases. The 
azimuthal velocity v<f, is significantly larger than the ra- 
dial velocity at r ~ 4i?*, but they become equal at larger 
distances. The sound speed c s is high inside the shock 
wave owing to the Joule heating J 2 /cr. Nevertheless, the 
outflow is supersonic in the region r ~ (4 — 7)ii*. Note, 
that above and below the equatorial plane the region of 
the supersonic outflow is larger (sec sonic surface line in 
the Figures 3 and 4). Figure 5b shows the variation of the 
velocities along the z— axis. 

Figure 6 shows the streamlines of the matter flow. Mat- 
ter free falls along the field lines going into the poles of the 
star. Some matter which flows close to the z— axis accretes 
to the surface of the star. However, matter more removed 
from the z— axis comes close to the star, it is deflected by 
the rotating magnetic field, and it then moves outward in 
the equatorial plane. 

Figure 7 shows the radial dependences of the different 
radial forces in the equatorial plane. One can see that 
for r > 3.5i?*, the centrifugal force becomes dominant 
in accelerating the matter outward. However, at larger 
distances, r > 5.5i?*, the pressure gradient force become 
larger and determines acceleration of matter. Thus, cen- 
trifugal and pressure gradient forces accelerate matter in 
the radial direction. Note, that in most of the region 
(r > 4.3i?*) the magnetic force is negative so that it op- 
poses the matter outflow. 

The results we have shown are for a relatively strong 
propeller. If the rotation rate is smaller, then the matter 
outflow become less intense. For < 0.16, no equato- 



rial outflow is observed and the shock wave becomes more 
nearly spherical as it is in the non-rotating case. On the 
other hand the shock wave becomes more a more flattened 
ellipsoid for larger w*. Matter rotates rigidly inside the 
Alfven radius ta while at large distances SI = v^/r de- 
creases. Figure 8 shows the radial dependences of the an- 
gular velocity in the equatorial plane for different values 
of w*. 

Figure 9 shows that only a small fraction of the Bondi 
accretion rate accretes to the surface of the star. Fig- 
ure 9a shows that the accretion rate to the star decreases 
as the angular velocity of the star increases, M/Mb ~ 
0.0075 w^ 10 . For a non-rotating star, M/M B w 0.07 in 
agreement with our earlier work (T99). Figure 9b shows 
that the accretion rate to the star decreases as the star's 
magnetic moment increases, M /Mb ~ fi ■ We have 
also done a number of simulation runs for different mag- 
netic diffusivities in the range f] m = 10~ 6 — 10~ 4,5 , and 
from this we conclude that M /Mb oc {rj m ) 0A . Accretion 
along the rotation axis of stars in the propeller regime 
was discussed by Nelson, Salpeter, and Wasserman (1993). 
Three-dimensional instabilities (e.g., Arons & Lea 1976) 
not included in the present two-dimensional simulations 
may act to increase the accretion rate to the the star's 
surfaces. 

Figure 10 shows the total angular momentum loss rate 
from the star as a function of w*. This was obtained by 
evaluating the integral 

L = -JdS-(pv p rv+-^^) , (9) 

over the surface of the numerical star. Here, dS is the out- 
ward pointing surface area element and the p— subscript 
indicates the poloidal component. The dominant contri- 
bution to equation (9) is the magnetic field term, and this 
term increases relative to the matter term as fi* increases. 
The dependence predicted by equation (3), L oc — f^ 2 ^ 13 
is shown by the dashed line. The fact that the observed 
dependence, L oc — fi^ 3 is steeper may be due to the fact 
that the magnetospheric radius is not much larger than 
the corotation radius. Also, we find that L oc — fi o s for 
a factor of four range of \i about our main case. From 
simulation runs with magnetic diffusivities in the range 
Vm = 10~ 6 — 10~ 4 - 5 , we conclude that L is approximately 
independent of r] m . 

5. ASTROPHYSICAL EXAMPLE 

Here, we give the conversions of the dimensionless vari- 
ables to physical values for the case (3 — 10 -7 , b = 10, 
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and fj m = 10~ 5 (equation 8). We consider accretion to a 
magnetized star with mass M — 1.4 M Q = 2.8 x 10 33 g. 
The density of the ambient interstellar matter is taken to 
be poo = 1.7 x 1CT 24 g/cm 3 (n^ = 1/ cm 3 ). We consider 
a higher than typical sound speed in the ISM, Coo = 100 
km/s, in order to have a smaller Bondi radius compatible 
with our simulations. The Bondi radius is 

R B = GM^/clc w 1.9 x 10 12 Clo 2 cm , (10) 

and the Bondi (1952) accretion rate is 

M B = 4 7 rA(GM JVS ) 2 Ax>A4 w 0.7xl0- 17 (n 1 /c 3 oo )M /yr , 

(11) 

where A = 0.625 for 7 = 1.4, ri\ = n^/l/cm 3 , and cioo = 
Coo/(100 km/s). Using the definitions (3 = Snpoo/B^, 
Poo = Pooclof-f, and the fact that b = 10, we obtain 
the magnetic field at the surface of the "numerical" star 
B, m 1.7 G. 

Recall that the size of the numerical star is R* = 
0.0044i? B w 0.82 x 10 10 cm in all simulations. Thus the 
magnetic moment of the numerical star is \x = B^r^/2 S3 
0.48 x 10 30 Gem 3 . Note that this is of the order of the 
magnetic moment of a neutron star with surface field 10 12 
G and radius 10 6 cm. 

The numerical star rotates with angular velocity f2* = 
oj^k*, where fix* = \/GM/i? 3 . For most of the fig- 
ures w* = 0.5, and this corresponds to a period of P* = 
2tt/0* w 685 s. 

Using the mass accretion rate of equation (11) and 
H30 — 1, the magnetospheric radius of the non-rotating 
star, from equation (1), is r m0 S3 6.6 x 10 10 cm £3 8i?*. The 
magnetospheric radius of the rotating star with uj* = 0.5, 
from equation (2), is r m Q S3 3.1 x 10 10 cm sa 3.8i?*. From 
Figure 8 note that the equatorial profile of the angular ro- 
tation rate Q(r, 0) for lj* = 0.5 is approximately flat out 
to rfi a t ~ 4.2i?* which is roughly the same as r m Q. As 

required, the corotation radius, r cor = i?*/^ 2 / 3 S3 1.6/?*, 
is appreciably less than r m n. For smaller values of uo* in 
Figure 8, r/; at decreases gradually. This dependence on uo* 
differs from that of equation (2) which may be due to the 
fact that r m Q is not very large compared with r cor . As lu* 
decreases, r cor increases rapidly. Thus the case a->* = 0.1, 
where r cor > rfi at , is not in the propeller regime. 

For the considered case of/3 = 10~ 7 , 60 = 10, f] m = 
10~ 5 , and uo* = 0.5, we have evaluated equation (9) numer- 



ically with the result that L -4.7 x 10 28 g(cm/s) 2 . The 
theoretical estimate of L from equation (3) for the men- 
tioned values of /i and r m o is L t h = 1-5 x 10 28 g(cm/s) 2 . 
A small, 70% reduction of r m n brings L t h into coincidence 
with the L from the simulations. 

6. CONCLUSIONS 

Axisymmctric magnetohydrodynamic simulations of 
Bondi accretion to a rotating magnetized star in the pro- 
peller regime of accretion have shown that: (1) A new 
regime of matter flow forms around a rotating star. Mat- 
ter falls down along the axis, but only a small fraction of 
the incoming matter accretes to the surface of the star. 
Most of the matter is expelled radially in the equatorial 
plane by the rotating magnetosphere of the star. A low- 
density torus forms in the equatorial region which rotates 
with velocity significantly larger than the radial velocity. 
Large scale vortices form above and below the equatorial 
plane. (2) The star is spun-down by the magnetic torque 
and to a lesser extent the matter torque at its surface. 
The rate of loss of angular momentum L is proportional 
to — f^ ,3 /i 0,8 , and it is approximately independent of r\ m . 
This dependence differs from the predicted dependence of 
equation (3) probably because the corotation radius is not 
much smaller than the magnetospheric radius r m . The ro- 
tational energy lost by the star goes into the directed and 
thermal energy of plasma. (3) The accretion rate to the 
star is much less than the Bondi accretion rate and de- 
creases as (a) the star's rotation rate increases (oc f2~ 10 ), 
(b) as the star's magnetic moment increases (oc /it -1 ' 7 ), and 
as the magnetic diffusivity decreases [oc (r] m ) 0A ]- (4) Be- 
cause the accretion rate to the star is less than the Bondi 
rate, a shock wave forms in our simulations and propa- 
gates outward. It has the shape of an ellipsoid flattened 
along the rotation axis of the star. 
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Fig. 7. — The figure shows the radial dependence of the radial forces acting to matter in the equatorial plane. The vertical 
scale is arbitrary. The centrifugal force is the dominant one which drives matter to the centrifugal equatorial wind. 
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Fig. 8. — Angular velocity of matter Vl — v^/r versus r in the equatorial region for different angular velocities of the star. 
The dashed line represents the Keplerian angular velocity VIk = \J GM/r 3 . 



Fig. 9. — The left-hand panel (a) shows the fraction of Bondi accretion rate reaching the star as a function of the star's 
angular velocity uo* = VI^/VLk*- For a non-rotating star, M/Mb = 0.07 or log(M/Ms) = —1.1. The right-hand panel (b) 
shows the dependence on the magnetic moment \i. The straight lines show the approximate power law dependences. 



Fig. 10.— Total angular momentum outflow from the star as a function of = Q*/Qk*- The straight line is a power 
law fit L oc — (fi*) 1 ' 3 . This is somewhat steeper than the predicted dependence of equation (3), L oc — (fi*) 12 ' 13 , which is 
shown by the dashed line in the figure. The vertical scale is arbitrary. 



This figure "RTTLfigl.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig2.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig3.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig4.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig5.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig6.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig7.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig8.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfig9.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



This figure "RTTLfiglO.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/astro-ph/0209548vl 



